Induced superconductivity in noncuprate layers of the Bi2Sr2CaCu208+5 
high-temperature superconductor: Modehng of scanning tunnehng spectra 
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We analyze how the coherence peaks observed in Scanning Tunneling Spectroscopy (STS) of 
cuprate high temperature superconductors are transferred from the cuprate layer to the oxide layers 
adjacent to the STS microscope tip. For this purpose, we have carried out a realistic multiband cal- 
culation for the superconducting state of Bi2Sr2CaCu208+4 (Bi2212) assuming a short range d-wave 
pairing interaction confined to the nearest-neighbor Cu d^2_y2 orbitals. The resulting anomalous 
matrix elements of the Green's function allow us to monitor how pairing is then induced not only 
within the cuprate bilayer but also within and across other layers and sites. The symmetry properties 
of the various anomalous matrix elements and the related selection rules are delineated. 



PACS numbers: 68.37.Ef 71.20.-b 74.50. -|-r 74.72.-h 



I. INTRODUCTION 

Scanning tunneling spectra (STS) of the cupratesi"— 
clearly show the presence of superconducting gaps and 
the associated coherence peaks. The 'leaking' of super- 
conductivity from the cuprate layers into the oxide lay- 
ers is a form of proximity effect^"—. A recent STS studj*^ 
finds that the magnitude of the superconducting gap or 
the pscudogap is not solely determined by the local dop- 
ing, but is also sensitive to the nearby nanoscale sur- 
roundings, raising the broader question as to how super- 
conductivity is transfered across various orbitals/sites in 
the cupratesJ^ In this connection, we have recently de- 
veloped a Green's function based methodology for carry- 
ing out realistic computation of scanning tunneling mi- 
croscopy/spectroscopy (STM/STS) spectra in the nor- 
mal as well as the superconducting state of complex ma- 
terials, where the nature of the tunneling process, i.e., the 
effect of the tunneling matrix element is properly taken 
into account. In our approach, all relevant orbitals in the 
material are included in a multi-band framework, and the 
tunneling current is computed directly for a specific tip 
position on the semi-infinite surface of the solid. An ap- 
plication to the case of overdoped Bi2212 was reported in 
Refs. |Tl| and[l2, where it was shown, for example, that 
the striking asymmetry of the STS spectrum between 
high positive and negative bias voltages arises from the 
way electronic states in the cuprate layer couple to the 
tip: With increasing negative bias voltage, new tunneling 
channels associated with d^^ and other orbitals begin to 
open up to yield the large tunneling current. The asym- 
metry of the tunneling current at high energies could thus 
be understood naturally within the conventional picture. 



without the need to invoke exotic mechanisms. Results 
of Refs. [Ill and [12 show clearly that the STS spectrum is 
modified strongly by matrix element effects as has been 
shown previously for angle-resolved photoemissior>i^, res- 
onant inelastic x-ray scatteringi^, and other highly re- 
solved spectroscopies fi^"—. 

The STM/STS modehng in Refs. M and [H is based 
on invoking the common assumption that the pairing in- 
teraction in cuprates is d-wave, involving nearest neigh- 
bor d^2_y2 orbitals of Cu atoms. Nevertheless, our com- 
puted STS spectrum reproduces, in accord with exper- 
imental observations, the superconducting gap and co- 
herence peaks at the position of the tip, even though the 
tip is not in direct contact with the cuprate layer. Our 
STS modeling scheme thus provides a natural basis for 
examining how the pairing interaction, which is limited 
to nearest-neighbor Cu dx2-y2 orbitals in our underlying 
Hamiltonian, gets transferred to other layers and sites. 

This article attempts to address these and related is- 
sues with the example of overdoped Bi2212. Central to 
our analysis is the concept of tunneling channels, which 
allows us to identify the contribution to the total tun- 
neling current from individual sites/orbitals in the semi- 
infinite solid. Moreover, we can distinguish between regu- 
lar and anomalous contributions to the tunneling signal, 
which arise from the corresponding matrix elements of 
the Nambu-Gorkov Green's function tensor. The anoma- 
lous channels are physically related to the formation and 
breaking up of Cooper pairs. In particular, matrix ele- 
ments of the anomalous Green's function can be used to 
monitor the contribution to the coherence peaks in the 
STS spectrum resulting from specific orbitals/sites in the 
material. In this way, we delineate how the pairing ampli- 
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minatcd by the BiO layer, which is followed by layers of 
SrO, CuOa, Ca, CuOz, SrO, and BiOJ^-^ The tunneling 
current is computed using a 2-\/2 x 2\/2 real space super- 
ccU consisting of 8 primitive surface cells with a total of 
120 atoms. The crystal structure is taken from Ref. [2l|. 
The STM tip is modeled as an s-orbital lying at the apex 
of the tip. The electron and hole orbitals included in the 
computations are: {s,Px,Py,Pz) for Bi, Ca and O; s for 
Sr; and {'is,d^z^-r^,dxy,dxz,dyz,dx^-y^) for Cu atoms. 
This yields 2 x 58 electron (spin up) and hole (spin down) 
orbitals in the primitive unit cell, or a total of 2 x 464 or- 
bitals in the simulation supercell. The Green's function 
is computed using 256 equally distributed k-points in the 
supercell which corresponds to 8 x 256 = 2048 k-points 
in a primitive cell. 

The multiband Hamiltonian in which superconductiv- 
ity is included by adding a pairing interaction term A 



FIG. 1: (Color online) Theoretical (green) STS spectrum nor- 
malized as shown in the figure is compared with the experi- 
mental (dashed blue) spectrum^ in optimally doped Bi2212. 
Regular (red) and anomalous (turquoise) parts of the com- 
puted spectrum are shown separately. All computations are 
based on Eq. ((6]). Coherence peaks arise only from the 
anomalous component of the Green's function. 



tude travels from the nearest-neighbor Cu-sites to other 
sites and orbitals within the cuprate plane as well as out- 
side to the second cuprate plane and to the BiO/SrO lay- 
ers. The symmetry properties of various matrix elements 
are analyzed and related selection rules are worked out. 

An outline of this article is as follows. The introduc- 
tory remarks are followed in Section II with an overview 
of the relevant methodological details of the underlying 
Hamiltonian and of our STS formalism. Section III dis- 
cusses proximity effects and is divided into several sub- 
sections, which address pairing amplitudes in various lay- 
ers. Section IV discusses selection rules and issues related 
to the symmetry of the gap through an analysis of the 
anomalous matrix elements. It is divided into consid- 
eration of on-site cases where the pairing orbitals lie at 
the same horizontal position, and to cases where these 
orbitals lie at other sites in the lattice. Finally, Section 
V presents a concluding discussion and a summary of 
our results. The added Appendix clarifies the symme- 
try properties of the regular matrix elements, which play 
an important role in the analyis of the symmetry of the 
anomalous matrix elements of the Green's function. 



II. DESCRIPTION OF THE MODEL 

The model underlying our analysis is the same as in 
Ref. [l2| to which we refer for details. An overview is nev- 
ertheless presented for completeness, and to introduce 
the various quantities needed for the present study. The 
Bi2212 sample is modeled as a slab of seven layers ter- 
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(1) 



with real-space creation (annihilation) operators cj^^. (or 
Ccct). Here a is a composite index identifying both the 
type of orbital and its site, and a is the spin index. Sa 
denotes the on-site energy of the a*'' orbital, and Vap is 
the hopping integral between the a and /3 orbitals. The 
hopping parameters are chosen to reproduce the LDA 
bands^-^ 

In the mean field approximation, the coupling between 
electrons and holes is of the form 



(2) 



ab 



Since the interaction U is not known, the standard prac- 
tice is to introduce a gap parameter, which gives the cor- 
rect gap width and symmetrjsSS. Specifically, we take 
A to be non-zero only between the dj.2-y2 orbitals of 
the nearest neighbor Cu atoms, and to possess a d- 



wave form, i.e., A 



d{d±x} 



lAI and A, 



d{d±y) 



where d denotes the (ij.2_y2 orbital at a chosen site, and 
dzt x/y the d^^-y^^ orbital of the neighboring Cu atom 
in x/y-directioni^ This form allows electrons of oppo- 
site spins to combine to produce superconducting pairs 
such that the resulting superconducting gap is zero along 
the nodal directions = i/cj,, and is maximum along 
the antinodal directions. The gap parameter value of 
|A| = AbmeV is chosen to model a typical experimental 
spectrum^ for our illustrative purposesi^i^ 

We discuss pairing between different orbitals in terms 
of the tensor (Nambu-Gorkov) Green's function Q (see, 
e.g., Ref. [M 



Ge F 



where Ge and Gt denote the electron and hole Green's 
function, respectively. 
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The following expressions for the pairing amplitudes in 
a tight-binding basis, which are derived in Ref. [T2| . are 
especially relevant for our analysis. 

(ca^c^t) = / de[l-2/(£)]pf^(e), (3) 

where the density matrix is 

pf^(e) = -i/m[F+5(£)]. 

Here, F^fjis) can be solved by using the tensor form of 
Dyson's equation for the retarded Green's function. Sim- 
ilarly, 

i^lAO- I ds[l-2fie)]pf^{e). (4) 

Eqs. ([3]) and (j4]) reveal the relationship between the 
anomalous part of the Green's function tensor and the 
pairing amplitudes between various sites. In particular, 
symmetry properties of Fa/s are seen to be related di- 
rectly to those of (c„|C^-|-). 

The tunneling spectrum is computed by using the 
Todorov-Pendry expression^^i^^ for the differential con- 
ductance a between orbitals of the tip (i, t') and the sam- 
ple (s, s'), which in our case can be written as 

^ = ^ = ^ E PtAEF)V,,,p,AEF+eV)V^l,. (5) 

tt' ss' 

Since electrons are not cigcnparticlcs in the presence 
of the pairing term, the density matrix can be rewritten 
by applying the tensor form of Dyson cquatioi>i^: 

= E(G-S"oG-,, -f F+^\F-A, (6) 

a 

where I]"q is the imaginary part of self-energy i^^'^^ The 
left-hand side of Eq. ([6]) is the ordinary density matrix for 
electrons, which is equivalent to the traditional Tersoff- 
Hamann approach^. However, as discussed in Ref. [T2| . 
our decomposition of the spectrum into tunneling chan- 
nels in Eq. (|6]) provides a powerful way to gain insight 
into the nature of the STS spectrum, especially in com- 
plex materials Note that the right side of Eq. (6) 
contains terms originating from the anomalous part of 
the Green's function. In Ref. [l^ we showed that coher- 
ence peaks appear only through the matrix elements of 
the anomalous Green's function. This role of the anoma- 
lous terms is demonstrated in Fig. [TJ where we see that 
the coherence peaks are absent in the partial spectrum 
resulting from the regular terms of the Green's function 
(red curve). 

III. INTERLAYER AND INTRALAYER 
PROXIMITY EFFECTS 

In this section, we analyze the induced pairing am- 
plitude (caicp-f) for a representative set of orbitals. It 



will be seen that despite the short range of the pairing 
interaction Aq,^, the anomalous Green's function, Fap, 
possesses a longer range. More specifically, we delineate 
induced pairing effects as follows: (i) Within a Cu02 bi- 
layer (Fig. HJa)); (ii) Intra-layer pairing in SrO and BiO 
layers (Fig. Hlb)); and (iii) Interlayer pairing between 
Cu02 and SrO/BiO layers (Fig. [H^c)). We discuss each 
of these cases in turn below with reference to Fig. [5] and 
Table H 



label 


orbital 


atom 


layer 







Cu (c) 


Cu02 (1st) 


1 




Cu (nn) 


Cu02 (2nd) 


2 




Cu (nn) 


Cu02 (2nd) 


3 


d,2 


Cu (c) 


CuOa (1st) 


4 


d,2 


Cu (nn) 


Cu02 (1st) 


5 


Pz 


0(c) 


SrO 


6 


Pz 


O (nn) 


SrO 


7 


Pz 


Bi (c) 


BiO 


8 


Pz 


Bi (nn) 


BiO 


9 


Px 


(b) 


Cu02 (1st) 



TABLE I: Shorthand notation used for the indices a and /3 
in Fig. [2] is defined. For each of the indices, varying from 
to 9, the table gives the atomic site [central (c), nearest 
neighbor (nn) and bonding (b)], the orbital and the layer 
involved. Order of the layers is: BiO, SrO, Cu02 (1st) and 
Cu02 (2nd), where BiO is the termination layer which lies 
closest to the STM tip. 



A. Pairing within Cu02 bilayer 

The most important anomalous matrix elements 
within the Cu02 bilayer are shown in Fig. [UJa). The 
red curve gives the contribution Fqi from d.j.2_y2 orbitals 
of two neighbouring Cu atoms in the x-direction, i.e., the 
matrix element between a spin up electron orbital at a 
Cu-site and a spin down hole orbital at the neighbour- 
ing Cu-site. This is the principal pairing matrix element 
since in our model is non-zero only between two such 
orbitals. 

The matrix element between the d,j.2_y2 orbitals of two 
Cu atoms at the same horizontal position within the 
Cu02 bilayer is zero by symmetry. However, the matrix 
element Fq2 between a Cu atom in the upper layer and 
each of the four neighbouring Cu atoms of the lower layer 
(and vice versa) is seen from Fig[2ja) to be substantial 
with an amplitude which is about l/4th of Fqi. This re- 
sult shows that pairing is not restricted to 6,^2 _y2 orbitals 
within a single Cu02 layer, i.e. it is not two-dimensional 
but extends vertically within the bilayer. 

The c?j,2 orbital of Cu also plays an important role. 
In fact, this orbital serves as a kind of gate for passing 
tunneling current from the cuprate layers to the SrO and 
BiO layers. Fig. [^Ja) shows that the amplitude F34 is 
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FIG. 2: (Color online) Imaginary part of the matrix elements of the anomalous Green's function Fa0 for various (q/3) pairs. 
[Recall that a and (3 are composite indices denoting both site and orbital.] Meaning of values of a and /3, which range from 
to 9, is explained in Table 1. For example, index refers to the d^2_y2 orbital on the central Cu atom in the cuprate plane 
nearest to the STM tip, and index 1 to the d^2_y2 orbital on nearest neighbor Cu atom in the same cuprate plane. The main 
matrix element between the two preceding orbitals, i.e., the (01) element, is shown by red lines for reference in all panels. Other 
matrix elements are shown scaled by factors ranging from 2-20 as indicated in the legends. Symmetries of the orbitals involved 
in various cases are shown schematically in the upper right hand side portions of the figures. Matrix elements compared with 
(01) (red) are: (a) (02) (green), (34) (dashed blue) and (99) (dotted black) for pairing within a Cu02 bilayer; (b) (56) (green) 
and (78) (dashed blue) for intralayer pairing in SrO and BiO layers; and (c) (03) (dotted black), (05) (green) and (07) (dashed 
blue) for pairing of Cu02 along the line connecting the central Cu and the surface Bi. 



about l/5th of Fqi and comparable to Fo2. There is 
also a smaller (about l/5th of rotationally invariant 
matrix element Fq^ (not shown in Fig. [2]) between the 
of a central Cu and the dj.2_y2 orbitals of the four 
neighbouring Cu atoms. At first sight this seems to break 
the d-wave symmetry, but we will show below that the 
combined symmetry of the orbitals involved remains d- 

42 



The role of 0-atoms in the cuprate layers can be de- 
lineated through the matrix elements Fgg and Fog. In 
Fig. [2][a) we show the onsite matrix element i^gg, which 
is about 1 / 4th of the Fqi term. We observe that in real 
space rotations of 7r/2 around the central Cu, Fgg changes 
it sign. A smaller contribution is found for Fgg (not 
shown in Fig. [2]). Its symmetry properties are consis- 
tent with the symmetry of the Zhang-Rice singlet, where 
a local orbital is constructed as a linear combination of 
the four oxygen atoms around the central Cu. The sym- 
metry analysis of Section IIVI below shows that both Fgg 
and -F09 are also consistent with the d-wave symmetry. 

Finally, we note that there is a substantial term, F03, 
between the (1^2 and dx2_y2 orbitals of the same Cu atom, 
which is perhaps surprising. Fig. [2{c) shows that this 
pairing amplitude is about half of Fqi. Since this is an 
onsite term, the d-wave symmetry again follows from the 
combined symmetry of the two orbitals, as discussed in 
Section IIVI below. 



B. Intralayer pairing within SrO and BiO layers 



In considering intralayer pairing in the SrO/BiO layers, 
we find that for Bi or apical O atoms, the most impor- 
tant non-zero anomalous matrix elements occur between 
Pz-orbitals of the central atom and its four neighbours, 
i.e. F56 and -FVs- These matrix elements possess the same 
d-wave symmetry as Fqi. While all matrix elements have 
the same energy dependence in Fig. ^h), Fqi is about 30 
times larger than F^q or Fyg. The coherence peaks lie at 
exactly the same energy in each layer, i.e., the gap width 
is the same in all layers. The scaling factor for the am- 
plitude seems to roughly follow the spectral weight of the 
orbitals. Hence, the pairing of electrons within the oxide 
layers seems to be a direct consequence of the tail of the 
Cu02 electron wave function within the various layers. 
This kind of pairing is in the spirit of the original idea of 
proximity effecl^r— where superconductivity is viewed as 
"leaking" from the superconducting part of the sample to 
the normal state material. Although the aformentioned 
orbitals are the most important ones, non-zero pairing is 
not restricted to just these orbitals. On the other hand, 
certain terms are strictly zero due to symmetry. In par- 
ticular, all the onsite Fs^, Fp^p^, Fp^p^, and Fp^p^ from 
Bi and 0(Sr) are zero, as are many Bi-O(Bi) and 0(Sr)- 
Sr terms. However, Fp^p of two neighbouring Bi's is 
non-zero as is Fp ,„ on the same Bi atom. 
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C. Interlayer pairing between Cu02 and SrO/BiO 
layers 

Pairing on BiO and SrO layers is not restricted to 
intralayer terms discussed above. The interlayer terms 
^05 and Fo7 between the p^-orbitals of the central Bi [or 
0(Sr)] and d,j.2_y2 of Cu right below these atoms is, in fact 
significant, while the anomalous term to the neighbour- 
ing Cu atoms is rather small. The existence of these ma- 
trix elements might be surprising, since the regular ma- 
trix elements are zero by symmetry between d^i_yi and 
the rotationally symmetric orbitals of the atoms above 
the central Cu (See Appendix R)) . But we will show in 
the following section that these matrix elements are not 
symmetry forbidden. From Fig. [2Ic), the scaling fac- 
tor between these elements and Fqi is of the order of 4. 
Notably, this interlayer pairing would appear as de- 
pendence in the gap-function. If we make a reflection of 
the slab with respect to the Ca plane lying between the 
two Cu02 layers, the corresponding anomalous matrix el- 
ements change their sign, indicating that this term has a 
node at fcj =0, and thus deviations from d-wave symme- 
try should be found for non-zero k^ - The fcz-dependence is 
also seen in the rather small terms between orbitals of 
the Bi atoms of the surface layer, and the nearest neigh- 
bour Bi atoms of the BiO layer half the primitive cell 
below the surface. 



IV. SYMMETRY AND SELECTION RULES 
FOR INDUCED PAIRING 

We now discuss the symmetry properties and the re- 
lated selection rules for the anomalous matrix elements 
of the Green's function in terms of the d-wave symmetry 
of the pairing matrix. 

A. Symmetry properties of the anomalous matrix 
elements 

Note first that the symmetry of the pair wave function 
depends on the relative motion of the pairing electrons, 
i. e., only on the relative coordinate Ri — Rj. The anal- 
ysis of the symmetry properties however becomes more 
transparent in fc-space. Accordingly, we transform the 
real space matrix elements Fiajp (e) into k-space as 

F„^(k, e) = ^(k|Oa)^^oa,^(e)(.7/3|k) 

' (7) 

j 

Here, we have set i?i = and (^^(k) is the orbital wave 
function in fc-space. The site indices i and j and the or- 
bital indices a and /3 are shown explicitly for all matrix 
elements. The summation is taken over the site index 
J. The orbital indices arc obviously not involved in the 



transformation. For simplicity, we will restrict the analy- 
sis below such that j is either on-site or one of the nearest 
neighbors of the central site. The generalization to far- 
ther out neighbors is straightforward. 

We need to take into account not only the phase differ- 
ence between the sites, but also the form of the tight bind- 
ing orbitals Lpaik). These orbitals have the same symme- 
try in real-space and fc-spacc. In particular, 

k 

= (kb.T) oc y 

^;(k) = (kb,) (X ^ 

/j2 _ ^2 (8) 
3fc2 _ ^2 

<^3^2_^2(k) = (k|d322_^2) CX . 

The symmetry of the matrix elements is now readily 
analyzed. We give two examples to illustrate the pro- 
cedure: (i) -Foi between the d^^-y^ orbitals of the cen- 
tral Cu and its neighbours; and (ii) i^03 between d^^ and 
dx2-y2 at the central site. In the case of _Foi the product 
of the orbital functions is even under rotations by 7r/2: 

^:(k)^„(k) = |(^„(k)p 

In fact, this applies to all cases where a — p. Summing 
over the four sites around the central Cu and applying 
the odd parity with respect to 7r/2 rotation of the real 
space matrix elements Foaja ^ ^oji we obtain 

FaaiK^) = 2|Foaja(e)| (cos (fc^fl) -cos(A:j,a)) |</3a(k)P, 

(9) 

which is obviously d-wave. This is easy to see for the 
dx2-y2 orbitals, but Eq. ^ leads to the same conclusion 
for any set of four neighbouring orbitals similar to the 
central one, as long as the real space element is odd under 
rotations by 7r/2. 

Turning to the case of F03, here the sum in Eq. ([7]) 
consists of a single term (central site), so that there is no 
site related phase factor. We only need to consider the 
product of orbitals: 

^*,2_,,2(k)^,2_,2(k) (X {kl - kl) (10) 

This is the only term through which angular dependence 
enters in Eq. ([7]). This again is d-wave, keeping in mind 
that only the in-plane symmetry is relevant. These con- 
siderations apply more generally to any case where the 
two orbitals involved lie at the same horizontal position 
with one of the orbitals being rotationally invariant and 
the other is d-wave. 



B. Selection rules for anomalous matrix elements 

Selection rules for the matrix elements Fiajp of the 
anomalous Green's function do not follow directly from 
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those for the corresponding regular matrix elements dis- 
cussed in Appendix A. For this purpose, we write Fi^jp 
asi^ 

F^cjpie) = -~G+^kyie)Ak^isGt,pi~e) (11) 

where the Einstein summation convention is implicit, and 
both the Green's functions on the right hand side of 
the equation are regular. The first is the renormalized 
Green's function for the superconducting state, while the 
second with superscript zero is the bare Green's function 
for the normal state. However, as shown in Appendix [Ki 
the symmetry properties of these two Green's functions 
are the same since both are regular. 

Equation ((TT|) highlights the central role of the pairing 
matrix Ak-yis in determining the symmetry properties of 
Fiajf}. However, summation over the intermediate states 
is cumbersome. Therefore, we convert Eq. (|lip to k- 
space first: 

F(k,£) = -G(k,e)AkG°*(k,-£) (12) 

where orbital indices are suppressed. Eq. makes 
it clear that F(k, e) possesses the d-wave symmetry of 
Ak since the regular Green's functions are rotationally 
invariant as shown in Appendix A. 

Converting F(k, e) to the real-space, yields 

^-,/3(e) = ^(*a|k)^'(k,e)(k|j/3), (13) 

k 

or 

i^oa.Me) =E'^'*''''^'"'''"'V:(k)^^(k)F(k,£), (14) 

k 

where we have fixed the first site index to Ri = 0, and 
made a separation into perpendicular (k^Zj) and parallel 
directions {k^Xj + kyi/j = fc|| ■ R\\j)- 

We first consider the case where orbitals a and f3 are 
at the same horizontal site, i.e., R\\j = 0. The necessary 
condition for the matrix element F^ajp to be non-zero is 
that 

^:(k)(^;5(k)F(k,£), (15) 

is rotationally invariant. Since i^(k, e) is d-wave, the 
product of wave functions on the right hand side of Eq. 
(|14p must have d-wave symmetry. For example, one of 
the orbitals could be d-wave symmetric and the other ro- 
tationally invariant. In particular, the anomalous matrix 
element between d^i_y2 and one of the set {s,Pz,d^2} 
satisfies this condition. 

Furthermore, if Zj =0, i.e., the orbitals are at the same 
site, the term (|15p must also have an even parity in the 
z-direction for a non-zero matrix element, and hence Pz 
would not be a possible pair with dx2-y2. However, in 
the case of Zj ^ 0, Pz is allowed, since the phase factor 
gtk^zj (^Qfjg j^Q^ have a well-defined parity. Note also that 
the anomalous matrix element between c?j.2_„2 orbitals 



at the same horizontal position (i?||j = 0) is necessarily 
zero, since the term is odd under rotations of 7r/2. 

We next consider matrix elements between orbitals at 
neighboring sites where, R\\j = xj = ±a, or R\\j = yj = 
±a. In this case, we will see that there will always be 
a non-zero matrix element with a properly symmetrized 
combination of neighboring wave functions, and the se- 
lection rules determine the correct choice of phase factors 
between sites. We discuss a particular case in detail as an 
exemplar. Specifically, let us compare sites Rj = (a, 0,c) 
and Ri = (0, a, c) and check whether or not the sign of 
the sum in Eq. ([T^ changes. For the first site we get 

Foajp = e'("^''+'=^^V:(fc., fcy)^^*(fcx, ky)F{-k, e) 

(16) 

and for the second site 

Foaw = J2 e'^'''''^''''''vl{kx,ky)^p{k,,ky)F{k, e). 

k 

(17) 

A rotation of 7r/2 is equivalent to the transformation 
ky kx and kx ^ ~ky Applying this to (|17p yields 

F^cifi = Y,<^'^''''"^'''"''^ai-ky^kx)<fip{-ky,kx)F{-k,e). 

(18) 

Thus the product ^p*f^{—ky,kx)^pp{—ky,kx) determines 
what happens under rotations of 7r/2 around the site 
1 = 0. There are two cases: (1) This product is equal to 
fai^x, ky)ip/3{kx,ky), so that these terms are invariant, 
and the total effect of rotation on F^ajp in Eq. fol- 
lows the d-wave symmetry of F(k, e); and (2) The prod- 
ucts of the orbitals in k-space have opposite sign, and the 
matrix element Foajp is invariant under in-plane rotation 
by 7r/2. In either case there will be pairing between the 
central orbital a and a properly symmetrized orbital 0, 
as defined in Eq. (jA6p of Appendix \K\ For case (1), 
an invariant 'fiaikx,ky)ipi3{kx,ky) linear combination of 
coefficients must be chosen with d-wave symmetry. For 
case 2, i.e., d-wave symmetric ip*^{kx,ky)ipp{kx,ky), the 
correct linear combination has all positive expansion co- 
efficients. Notably, for a — (5, the product of orbitals 
is invariant, so that any pair involving the same orbitals 
at neighboring sites must involve a linear combination of 
neighbors which is odd in rotations by 7r/2. For exam- 
ple, in the anomalous matrix element i^jg between the 
Pz-orbitals of two Bi neighbors, the coefficient in the a;- 
direction has an opposite sign to that in the j/-direction. 

V. DISCUSSION AND CONCLUSIONS 

We emphasize that the logic of symmetry rules for the 
anomalous matrix elements is more complicated than 
that of the regular matrix elements. In particular, 
the nonvanishing tunneling channels can be determined 
through group theoretic consider ation o ^ ^ i ^ ^ . For exam- 
ple, since the rotational symmetry of the Pz orbitals of Bi 
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and apical oxygen atoms differs from that of the d^^-y^ 
orbital of the Cu at the same horizontal position, the 
corresponding off-diagonal term of the regular Green's 
function vanishes, inhibiting the corresponding tunnel- 
ing channel. In contrast, coupling between electron and 
hole degrees of freedom via the gap matrix Aq^ leads to 
less obvious symmetry rules for the anomalous matrix el- 
ements: Now the quasiparticlcs are linear combinations 
of spin up electrons and spin down holes, and there is 
no simple rule for selecting the orbitals contributing to a 
chosen quasiparticle state. Hence, the Pz or (1^2 orbitals 
of Bi, O or Cu atoms may couple to a c?^2_y2 orbital of 
a Cu atom at the same horizontal position, and the pos- 
sibility of this coupling must be checked by considering 
the tensor form of Dyson's equation, as written out in Eq. 
p2p , together with the transformation into tight-binding 
basis of Eq. (fM)) . 

In summary, we have presented a comprehensive study 
of anomalous matrix elements of the Green's function de- 
rived from a realistic multiband model of Bi2212. The 
imaginary parts of these matrix elements describe the 
contributions of different orbitals to the coherence peaks 
involving the formation and breaking up of Cooper pairs. 
Although the pairing interaction is modeled by a local d- 
wave term in the Hamiltonian connecting only the dx^-y-^ 
orbitals of neighbouring Cu atoms, the anomalous matrix 
elements display a longer range with induced supercon- 
ductivity appearing at other sites/orbitals, including the 
second cuprate layer and the BiO/SrO overlayers. Our 
analysis delineates the precise routes through which the 
induced superconductivity in a complex cuprate system 
is transferred between various orbitals and sites. 
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Appendix A: Symmetry properties of regular matrix 
elements of the Green's function. 

This appendix delineates the symmetry properties of 
the regular matrix elements of the normal and super- 
conducting (SC) state Green's functions G°(e) and G[e), 
respectively, which were seen in connection with Eq. 
above to be important for understanding the nature of 
anomalous matrix elements. Taking the origin at the po- 
sition of the Ot/i atom, Gqi^^^(£) can be written as 

(Al) 

k 



where 

G°(k,£)= ^ ^ ^ (A2) 

Since the Hamiltonian is invariant under rotations of 7r/2, 
the dispersion and G'°(k, e) are also invariant. This is 
true as well for the SC state regular Green's function 
since the self-energy in Eq. \K2\ is augmented by an 
additional term Y,^^^ = AfcG5|(k, £)A|,, which is rota- 
tionally invariant^. Because G° and G possess the same 
symmetry properties, in the following, we only consider 
the symmetry properties of G°(k, e). 

Consider first the case where Rj = (0,0, c). Then, 
Goajp 7^ only if (^* (k).^^(k) in Eq. (|1T|) is invariant 
under the in-plane operations of the symmetry group of 
the Hamiltonian. For example, a Pz orbital can have 
non-zero matrix elements with s, pz or dz^ of an atom 
at the same horizontal position. But the matrix element 
between pz and d^^^y-^ of atoms at the same horizontal 
position is zero. For c = 0, the matrix element is non-zero 
only if the orbitals are similar. 

We next consider the case where there are four atoms 
around a central atom at the distance of the horizon- 
tal lattice constant a : Ri = (a, 0,c), i?2 = (0, a,c), 
i?3 = (— a, 0, c), and i?4 = (0, — a,c). Changing the in- 
dices according tol^2^-3^4^-l corresponds 
to rotations by 7r/2 in real space. The transformation 
ky — >■ kx and — > — fcj, represents the same rotation in 
k-space. Now the phase factor e**''^J has a fundamental 
effect on the symmetry behavior. Let us compare cases 
Ri = (a, 0, c) and R2 = (0, a, c) and check whether the 
sign of the sum changes. In the first instance we get 

G'Li/3(e) = Ee'^'^"^'^'V:(fc.,A:,)^0(fc,,fc,)GO(k,e) 

(A3) 

while the second case gives 

Gl^2pi.e) = 5]e'('=="^+'=--V:(fc.,A:,)^^(fc„A:,)G°(k,e) 

(A4) 

Applying the transformation ky — > k^ and k^ ^ —ky to 
(|Xt)) yields 

Gl^e) = e'^'-^'^'^'^'pli-ky, kx)M-ky, kx)G'{K e). 

(A5) 

Thus, it is the product ip'^{—ky,kx)ipp{—ky,kx) that 
determines what happens under rotations of 7r/2. If this 
is equal to ip*^{kx, ky)ipj3{kx, ky), the matrix element does 
not change sign under rotations, otherwise Gq^^^ changes 
sign under in-plane rotation of tt/2. In particular, for 
a = /3, this term is invariant, but for a — pz or dz^ and 
/3 = dx2_y2, there is a change of sign under rotation. 

An equivalent approach is to consider a linear combi- 
nation of orbitals jP 

4 



There is a non-zero regular matrix element Gg^ ^ only if 
belongs to the same representation of the symmetry 
group of the Hamiltonian as orbital |0a). The transfor- 
mation of the expansion coefScients Cj directly follows 
from the transformation of G^^^^. For example, it is ob- 
vious that 

4 

Hence, if a = /3 = dx^-y^ and j = 1...4 are defined as 
above, Cj must be a constant in order to keep the full 
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symmetry of the group of the Hamiltonian, leading to 

Gla,4> cie*'=^'=[cos (fc^a) -f cos [kya)]{kl - k^f 

If, however, a = d^-i and (3 = o?^2_,y2, one must have 
C2 = C4 = — Ci = — C3, and then 

GL,^ cx cie^'^-^icos (fc,a) - cos (feya)](fc2 - kl){3kf^ - fc^), 

which requires change of sign of Cj 's under rotations of 
tt/2 in order to obtain an invariant matrix element. 
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